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^ ! Abstract 

> . 

2 We study in detail the application of renormalisation theory to models of cluster aggregation and frag- 
^ mentation of relevance to nucleation and growth processes. We investigate the Becker-Doring equations, 
'"O originally formulated to describe and analyse non-equilibrium phase transitions, and more recently gen- 
^ \ eralised to describe a wide range of physicochemical problems. In the present paper we analyse how 
q ■ the systematic coarse-graining renormalisation of the Becker-Doring system of equations affects the ag- 
L " ~" '. gregation and fragmentation rate coefficients. We consider the case of power-law size-dependent cluster 
t-h ! rate coefficients which we show lead to only three classes of system that require analysis: coagulation- 
^ dominated systems, fragmentation-dominated systems and those where coagulation and fragmentation 
■ are exactly balanced. We analyse the late-time asymptotics associated with each class. 
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1 Introduction 



The purpose of this paper is to develop renormalisation methods for nucleation and growth processes. Renor- 
malisation theory has been applied to various problems in theoretical physics, being widely used in high energy 
physics ^3] and in the analysis of equilibrium phase transitions in statistical physics 0. It has also been 



used in a limited way in asymptotic analysis by Goldenfeld |yj . However, it has not previously been used in 



the asymptotic analysis of non-equilibrium phase transitions of the kind to be discussed in the present paper, 
mainly due to the complex nature of these models which often involve nonlinear and non-local features. 

In statistical mechanics, the basic idea underlying the concept of renormalisation is the systematic filtering 
out of unwanted degrees of freedom in the transition from the microscopic to the macroscopic description 
of phenomena. In the case of equilibrium phase transitions, near a critical point the system looks the 
same on all length scales and this physical insight is translated into mathematical language in terms of 
'the renormalisation group' (RG), a set of scale-invariant transformations which leave the essential physical 
properties of the system unchanged. Admittedly, the RG terminology is somewhat inappropriate since the 
symmetry of such transformations is that of a semi-group (information is lost on coarse-graining), and the 
methodology comprises many different ideas and distinct methods, rather than being a formal monolithic 
edifice. Wilson's work on the renormalisation group spawned a vast effort in the analysis of equilibrium 
critical phenomena and quantum field theory |38|, |I| . Of much more recent interest has been the application 
of the renormalisation methodology to non-equilibrium phenomena. The complexity of far from equilibrium 
dynamics is such that, while the physical motivation behind the RG programme of coarse-graining microscopic 
models still seems appropriate to obtain the macroscopic properties, it must be shown in each specific case 
that there may be a suitable scale-invariance ("renormalisation group") underlying the physics, which can 
reveal or explain universality in behaviour. 

Of central interest in the present paper is the asymptotic late time macroscopic behaviour of complex 
dynamical systems. Previously, Bricmont & Kupiainen have combined asymptotic methods with ideas from 
renormalisation theory to study problems in diffusive processes and special limiting solutions of nonlinear 
parabolic equations [Q, |], [| while Woodruff has cast multiple timescale problems in the form of renormalisa- 



tion theory [|39] , |42| , [40], |41] . This allowed the equations governing larger-scale phenomena to be derived from 
a more general theory and to be separated out from the small scale dynamics. In his approach, Woodruff 
separates timescales in the asymptotic limit e — > by writing the time t as t — t/e + 1 where t is the fast 
timescale and t is the slow timescale. The methods we develop in this paper share some of these features. 

We apply underlying concepts from renormalisation theory to study a simple model of nucleation, namely 
the Becker-Doring equations. These were originally formulated to study first-order phase transitions, and have 
been a central pillar in classical nucleation theory. They describe the stepwise growth and fragmentation of 
clusters in terms of the rates of the individual processes wherein monomer particles join or leave each cluster. 
Despite the age of this model, we have only very recently applied generalisations of the Becker-Doring 
equations to a wide range of physicochemical processes from those involving surfactant self-assembly |14|, [3j| 
through RNA-polymer- formation JT^L [T3|, [3^] to cement-setting jnS, |3"I|| . In these studies coarse-graining 
procedures are derived and used to reduce large systems of equations down to progressively lower-dimensional- 
"mesoscopic" -dynamical systems capable of theoretical analysis using standard techniques from the modern 
theory of differential equations. The coarse-grained contraction procedure, which is summarised in Section 2, 
is analogous to other renormalisation methods used in statistical physics, and this feature is exploited within 
this paper. 

The Becker-Doring equations have recently been subjected to more conventional analysis using matched 
asymptotic expansions [56 , ff7|| . Rather than use the coarse-graining approximation, which emphasises the 



discrete nature of the equations, that analysis concentrated on the large-time limit where continuum ap- 
proximations become valid [p6 |. The key results from this analysis are summarised here in Section 3. The 
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large-time limit of the coarse-grained equations was also studied there, and revealed the necessity of combin- 
ing a change of timescale with the coarse-graining process in order to preserve the accuracy of the dynamics of 
the system; in the present analysis, this translates into the need to use a dynamical renormalisation scheme. 

Asymptotic methods have been applied directly to coagulation-fragmentation problems before, for exam- 
ple van Dongen & Ernst Hl8 , [T9[1 and Davies et al. [|17[] have elucidated the self-similarity of the pure aggrega- 
tion form of Smoluchowski's equations. These equations differ from the Becker-Doring equations in allowing 
clusters arbitrary size to coalesce together whereas the Becker-Doring equations only allow cluster-monomer 
interactions, but allow both aggregation and fragmentation. The Smoluchowski coagulation-fragmentation 
problem is much more complex but has been analysed by van Dongen & Ernst |2(| and Carr & da Costa 
H [5| |TIJ. The scaling behaviour of the Becker-Doring equations has been studied by several authors, for 
example, Brilliantov & Krapivsky || and Blackman Sz Marshall [||. Coarse-grained versions of the Becker- 
Doring equations have been analysed using matched asymptotic methods |[30|| , [ |33|j . In addition many authors 



have considered hybrid systems which have combined features of the Smoluchowski equations and the Becker- 
Doring equations: this approach started with Samsel & Perelson who studied a system in which two 
clusters of arbitrary size can coalesce, but with a fragmentation step which is strictly Becker-Doring in na- 
ture, allowing only monomers to dissociate from clusters. Such models have been investigated further by 
Krapivsky & Redner |25| and Hendricks & Ernst P^| . Krapivsky |24| considered an alternative modification 
of the Becker-Doring equations, again concentrating analysis on the large-time scaling laws, while Peliti [26] 
looked at the scaling behaviour of a very simple model of pure aggregation in various dimensions in order to 
discover the manner in which the scaling laws depend on spatial dimension. 

Our methods fit within the spirit of Woodruff's approach if we write the microscopic aggregation number 
r of a cluster as r = (A — l)n + 1 + k where n is of mesoscopic size and A; is a microscopic correction; we 
then aim to determine the problem on the mesoscale in a form which does not require us to simultaneously 
solve the microscopic problem. Thus microscopic detail is filtered out, but we are able to construct a simpler 
model which remains valid on larger (mesoscopic and macroscopic) scales. Velazquez has recently used a 
renormalisation technique in an attempt to draw together the theories of Lifshitz-Slyozov coarsening and 
nucleation as modelled by the Becker-Doring equations [29|. His renormalisation group differs from ours in 



several ways, most egregiously in that it separates the monomer concentration from the large-cluster sizes, 
whereas we retain the monomer concentration as an integral part of the mesoscale model. 

Preliminary results of our work were reported in an earlier publication |16|| . Here we provide the details 
omitted from that communication. We also discuss how the coarse-graining procedure affects the large-time 
kinetic behaviour of the systems under consideration. Various types of system are introduced and analysed, 
including full and coarse-grained systems of Becker-Doring equations (Sections |2| and |3[), the underlying renor- 
malisation structure described (Section £|), and the large time asymptotics associated with the fixed points 
of the renormalisation mapping are analysed (Section We consider noise- free Becker-Doring systems, by 
which we mean those with simple analytical formulae for their rate coefficients; in a separate paper |34| we 
present our analysis of "noisy" systems in which the rate constants are perturbed by a random amount from 
the analytical expressions used here. We conclude with a discussion of our results (Section^). 



2 The Becker-Doring cluster equations 

In this section, we give a basic outline of the Becker-Doring system of equations and their properties. We 
consider the case for which the monomer concentration (ci) is held constant; thus the Becker-Doring equations 
we are concerned with here are 

c r = J r -\ — J r , (r > 2), J r = a r c r C\ — 6 r+ iC r+ i, (2-1) 
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and we leave the constant mass formulation (in which the monomer concentration may vary) for future study 
In this system of equations, c r (t) represents the concentration of clusters containing r monomers. Clusters 
of any particular size can be formed by two processes: either by the next smallest cluster size coalescing 
with a monomer, or by the next largest size losing a monomer. Only such monomer-cluster interactions are 
permitted in the Becker-Doring model of nucleation; cluster-cluster interactions are ignored. The assumption 
of a constant monomer concentration is made in situations where the so-called "pool chemical approximation" 
is valid, namely where there is a large source of monomer species entering into the system at a rate which 
maintains the monomer concentration essentially fixed and independent of time. 

There are certain special properties of the Becker-Doring system which are not immediately apparent 
from the equations ( |2.1| ), and which must be preserved under the coarse-graining rescaling if the method is to 
produce physically relevant results, through approximating the original system in a faithful manner. These 
are: 

(i) There exists a partition function, Q r , satisfying a r Q r = b r+ iQ r+ i, and hence an equilibrium solution 

c r = Q T c\. The quantity Q r is related to the chemical potential of a cluster of size r. If we denote 
Boltzmann's constant by k, the temperature of the system by T and use the reference state in which 
[if = 0, then the chemical potential of a cluster of size r (denoted /i r ) is given by [i r = fif + kT\ogc r . 
At equilibrium, this must satisfy /i r = 77/1, thus fif = —kT\ogQ r . 

(ii) There exists a Lyapunov function, V = J2^Li c r (}og(c r /Q r c[) — 1), which is related to the Helmholtz free 

energy of the system by F = kTVV, where V is the volume of the system. 

(iii) There exists a set of identities (or "weak form"), g r c r = g\J\ + [flv+i — g r ] Jr- 

r=2 r=l 

oo 

(iv) The total density of the system is g = ^rc r , however this is not conserved since matter can be added 

r=l 

to or removed from the system by he monomer concentration being held constant. 

For certain choices of rate coefficients a r ,b r and certain monomer concentrations ci, the equilibrium 
solution will not decay to zero in the limit r —>■ 00. In these cases, an alternative steady-state solution will 
be approached in the large-time limit. This solution is given by a constant nonzero flux through the system, 
that is J r = J independent of r. This condition yields the family of solutions 

Cr = QrC[ (l ~ j£ ^ fc+1 j , (2-2) 

\ fc=i a fcVfcCi / 

which contains the equilibrium solution as the special case J = 0. The steady-state flux J is determined by 
requiring the concentrations c r to decay to zero in the large-r limit, giving 



2.1 Other formulations 

An alternative formulation of the Becker-Doring equations has been proposed by Penrose & Lebowitz |27|j , 
in which the monomer concentration is allowed to vary, whilst the total mass of material in the system (g) is 
held constant. The analysis of such a system is in general more complicated since the monomer concentration 
is unknown and causes the system to become both nonlinear and non-local (in r). 
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An underlying model in the present work is a constant mass system in which monomer is added to or 
removed from the system by way of a precursor chemical P which undergoes a reversible reaction to create 
monomer C\. This is modelled by the simple reaction P ^ C% with forward rate kf{p,C\) and reverse rate 
kb{p, Ci). Denoting the concentration of P by p(t) we then have the system 

p = k b (p, ci)c! - k f (p, d)p, (2.4) 

oo 

ci = kf(p, c{)p - h{p, ci)ci - Ji - ^2 J r , (2.5) 



r=l 



in addition to ( |2.1| ). This system has a conserved quantity g tot = p + Y^°=i rc r-, however note that the 
amount of mass in the Becker-Doring part of the system ({cv}^) is not constant £?bd(£) = Y^°=i rc r- Such 
an extended system of the Becker-Doring equations has been used previously to model generalised nucle- 
ation phenomena |3l|]. In certain cases analytical progress in the study of such systems can be made using 
asymptotic techniques, such as in where the constant mass formulation of the Becker-Doring equations 



are investigated, and in |3l[ where a generalised constant mass system with inhibition as well as a precursor 
chemical is analysed. 

In studying just the Becker-Doring part of this extended system ( |2.1| ), we are assuming that the precursor 
chemical decays in such a way that the monomer concentration remains constant, that is c% =0. This requires 

oo 

P = - h - £ J r- ( 2 - 6 ) 

r=l 

Solutions which we describe as equilibrium solutions satisfy J r = for all r > 1, thus correspond to p — 
ki,ci/kf, and for this specific concentration of precursor chemical the additional reaction which adds or 
removes matter from the Becker-Doring part of the system will also be in equilibrium. Thus the equilibrium 
solution of our extended system is given by c r = Q r c\ with c\ = kjp/k b , and at equilibrium J r = for all 
r > 1. However for solutions we describe as steady-states, the fluxes J r = J ^ where J is independent ot 
size (r), so there is a net flux of matter into the system (assuming J > 0, consequently a removal of mass if 
J < 0), and in these cases we assume that there is sufficient precursor chemical to sustain the steady-state. 
In this case although c r = for all r, p ^ so the precursor chemical is not in equilibrium. Such steady-state 
solutions can be characterised by the flux associated with them. 



3 Coarse-graining procedure 

Following a general coarse-grained contraction with constant "mesh" size A in aggregation number (so that 
we only retain the aggregation numbers r = A n = (n — 1)A + 1), the kinetic equations reduce to 

in = L n _i — L n , (r > 2), L n = a n x n x^ — (3 n+ ix n+ i, (3-7) 

a n = Ta An a An+1 . . . a An+1 -i, /3 n+ i = Tb An+1 b An+2 . . . b Ari+1 , (3.8) 

where the retained coarse-grained cluster concentrations are relabelled with X\ := Ci the monomer 

concentration, which is not involved in the coarse-graining since it has a special role in the Becker-Doring 
theory fli~4 , 16| . The parameters a n , j3 n are the coarse-grained aggregation and fragmentation rates, now 



representing the addition or removal of A monomers to or from a cluster (rather than just a single monomer 
which occurs in the full Becker-Doring system). This flux of matter is denoted by L n . In effect, the system 
models the aggregation-fragmentation processes 

(X + 1)X 1 ^X 2 , X n + \X^X n+l (n>2). (3.9) 
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for coarse-grained clusters of size X r . The concentration x n (t) is representative of the concentrations c r for 
cluster sizes (A n _i + 1) < r < A n . The factor T represents a change of timescale. This factor was omitted 



in our original formulation of the coarse-graining procedure |14], |32| ; its inclusion ensures that the large time 
asymptotic behaviour of the reduced system coincides exactly with the original fine-grained system in the 
case of size-independent aggregation and fragmentation rates (a r — a, b r — b). The large-time asymptotics of 
the constant-monomer-concentration Becker-Doring system with constant coefficients were found by Wattis 
and King using conventional asymptotic expansion methods ||36|| . The leading-order solution found there 
falls into two cases: for ac\ < b the equilibrium solution c r = (a/6) r_1 c^ is approached; whilst for ac\ > b 
the steady-state solution c r = c\ is approached. In both cases the asymptotic solution is realised by a front 
moving forward into larger r-space ("aggregation space") at speed s = \ac\ — b\, ahead of which (r > s(t)) 
the concentrations are all zero, and behind which (r < s(t)), the concentrations are either at a steady-state 
or in thermodynamic equilibrium. In order for the coarse-grained contraction to preserve the correct speed 
of propagation, the timescale T should take the value 



T 



(aci — b) 
\(a x c\ - b x ) ' 



(3.10) 



This temporal rescaling also ensures that the steady-state fluxes in the full and reduced systems coincide. In 
the terminology of renormalisation theory, this makes the coarse-graining procedure a dynamical renormali- 
sation transformation. 

The four properties listed in Section |2] above for the full Becker-Doring system of equations are shared by 
the contracted system of eqns (|3.7f ): 



1 7, 



Qh.n X \ n ' 1 



oo oo 

9nX n = 9\Ll + J2(dn+1 ~ 9n)L n , 

n=2 



n=l 



V = x n ( log 



n=l 



X >■ 



QAn X l 

oo 

[{n- l)A + l]z n , 



n=l 



(3.11) 



where a n Q An = (3 n +iQk n+l - 

Note that it is also possible to allow the mesh size to vary with aggregation number [pi], that is putting 
A = A n transforms the system with constant rates (a r = a, b r = b) into one with non-constant rates (since 
a r = a Xr will be a nontrivial function of r, and similarly for (3 r+ i). For this reason we shall put to one side 
such generalised meshes in the present paper. 



4 Renormalisation structure 

We now proceed to consider the renormalisation of Becker-Doring models in which the cluster rate coefficients 
are of power-law form. As discussed below, this dependence is of immediate relevance to the description of 
surface-limited nucleation and growth processes in physicochemical systems. 

4.1 The Becker-Doring system with power-law coefficients 

In the case of Becker-Doring systems with power-law rate coefficients, the rate coefficients for aggregation 
and fragmentation are respectively 

a r = ar p , b r+ i = br p , (4.1) 

so that the parameter 6 = ac\jb which arises in the const ant- monomer case remains useful for classifying 
dynamical behaviour. The parameter p determines the variability of rate with cluster size, with p > 
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implying that large cluster sizes have larger aggregation and fragmentation rates, and p < giving rates 
which decrease with increasing cluster size. The latter case is of rather less physical relevance, but we 
still study it here for the sake of completeness. Typical values for p are p = 0, |, |, |, 1 for the examples 
of linear chain polymerisation, coagulation kinetics in two space dimensions, diffusion-limited coagulation 
in three dimensions, surface-limited coagulation in three dimensions, and branched chain polymerisation, 
respectively. Since a cluster's volume scales with aggregation number r, if we assume that clusters are 
spherical then their surface area scales with r 2//3 and their diameter with r 1//3 , accounting for the presence of 
these exponents. More general exponents can be manifest in other physical situations |TT . 



The partition function Q r is defined by Q r = {a/b) r 1 as in the p = case considered earlier [ 36| . The 
governing equations are then 

c r = J r -i — J r (r > 2), J r = r v {ac\c r — bc r+ i), (r > 1). (4.2) 

For 9 < 1 the system approaches the steady-state solution (which we shall also refer to as the equilibrium 
solution, even though strictly speaking it is not, owing to the boundary condition of constant monomer 
concentration imposed by the pool chemical approximation) given by solving J r = 0, that is 

c r = 9 r - 1 c 1 . (4.3) 

Note that this solution is independent of p, although the way in which the steady-state solution is approached 
depends on p. 

For 9 > 1 this solution gives diverging concentrations for large r and is hence unphysical. In this case the 
system approaches one of a family of time-dependent solutions to eqns ( |4.2[ ), in which all fluxes are equal; 
J r = J for all r then implies 

<v=^(i-/|:^). (4.4) 

For 9 > 1 the sum is convergent in the limit r — ► oo, and so the flux which gives the least singular behaviour 
in this limit is 

/oo 
k- p 9- k . (4.5) 
k=l 

The steady-state solution valid for all r is 

T OO 

c r = ^Y,k- p 9 r - k -\ (4.6) 

4.2 Repeated application of contraction procedure: renormalisation 

We have already shown how the full Becker-Doring system of equations can be reduced by a coarse-graining 



process to the subsystem of equations (|3T7|) -( |378|) . We now show that these equations can be further reduced, 



by applying the contraction procedure to ( |3.7| )-( |3T8D with mesh size fi, defining the monomer concentration in 



the new variables by Z\ = Xi(= C\). The new cluster concentrations z p are representative of the concentrations 
x n for (M p _! + 1) < p < M p where {M p } ! ^ =1 is the grid of retained n values of cluster aggregation numbers, 
determined by M p = (p — l)p, + 1. Thus we find 



z p — I p -i — I p , (p > 2), I p — ApZpZ^ — B p+ iz p+ i, (4-7) 

where I p is the flux of material from concentration z p to z p+ x, and has aggregation and fragmentation rates 
determined by 

A p = TaMpO-Mp+i ■ ■ ■ OiM p+1 -i, B p+ i = T/3 Mp+1 p Mp+2 . . . Pm p+1 - (4.8) 
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We have denoted the change of timescale by T. In the case of the model with constant coefficients (a n = a, 
(3 n = (3 for all n), in order for the timescales of the approach to steady-state or equilibrium to be correctly 
reproduced, the temporal rescaling factor T should be chosen to be 



~ = (ax x - (3) = (ad - b) 



by analogy with (|3.10|) . The rate coefficients satisfy v4 p (5( p _i)a m +i = B p+ iQ p x^ + \. The system satisfies the 
four properties 



p 



P =i 



,Q{p-l)X(j.+l. 



J2dpZp = 9ih + J2 (9p+i ~ 9p) 1 p Q = X][(p-l)V + l' 

p=l r=l p=l 



(4.10) 



So we see that the repeated contraction is identical to a single contraction with mesh parameter 7 = A/i. 
The system of equations ( |4.7| ) could have been derived by a single coarse-grained contraction, with mesh 

2p = c rr , zi = c x , r p = (p- 1)7 + 1. (4.11) 

Note also that T p = Am p = M\ p , so that the coarse-graining procedure is commutative. 

The combined effect on the rate coefficients and timescale can be verified by noting that in the case of 
constant coefficients, A = Ta^ 1 = TT^a x ^ = Ta Xtl , where 

- aci — b , . 

T = \ , 4.12 

which is what one would expect if A were replaced by A/i in equation (|3.10| ). 

At the heart of renormalisation theory is the repeated application of a coarse-graining procedure. The 
analysis above shows that: (i) a contracted version of the Becker-Doring system can itself be contracted, 
which is crucial if we are to use ideas from renormalisation theory; and (ii) the effect of coarse-graining the 
Becker-Doring equations with a mesh size A and then coarse-graining the resultant equations with a mesh 
size \x is identical to a single application of the coarse-graining process with mesh 7 = A/i. Thus the limit 
of repeated applications of the renormalisation coarse-graining process corresponds to the A ^> 1 limit of 
eqns (|3.7[). This second property makes the ensuing analysis possible, since we do not have to handle systems 
of equations more complicated than the Becker-Doring system itself; that is, the structure of the equations is 
invariant under this transformation. Note that as one may have expected a priori, the coarse-graining does 
not introduce cluster-dimer or more general cluster-cluster interactions; rather the theory is generalised so 
that clusters interact with several monomers in moving from one coarse-grained set of cluster sizes to the 
next. 



4.3 Contraction of power law coefficients 

If the rate coefficients in the original formulation in eqns ( |2.1| ) are determined by simple power laws, namely 
a r = ar p and b r+ i = br p , then the coefficients in the reduced model are 

a n = a x {[(r - 1)A + l][{n - 1)A + 2] . . . [n\]} p . (4.13) 



S 



Thus 



A 

loga n = A log a + p log(nA 
= Aloga + pA [log(nA) 

For large n this asymptotes to loga n ~ A log a + pA log(An), so for simplicity we shall take a n = (aX p n p ) x , 
which differs slightly at small values of n. The backward rate coefficient is then (3 n+ i = (b\ p n p ) x . Since 
contracting a constant density Becker-Doring system leaves the kinetic part of the equations invariant and 
only affects the coefficients in the "constitutive equation" for the flux, the steady-state solutions and large- 
time asymptotics for the contracted systems of equations can be found directly from the above theory. 

The new system has its own 9 parameter determining the balance between aggregation and fragmentation 
rates in the system which, for the moment, we shall call 9 = a r x x /(3 r+ i = 9 X , so contraction of the system 
maps 9 to 9 X . The parameter 9 thus plays an important role in the renormalisation procedure, the fixed 
points of this mapping corresponding to 9 = 0, 1, oo; hence such systems are of special interest to us. Also 
the contraction maps coefficients with exponent p to those with exponent pX. Thus following a contraction, 
there are only three limits to consider: small p (namely p — 0) and large p (positive and negative). 

The effect of coarse-graining a Becker-Doring system is to modify the rate coefficients, by the map 9 i— > 9 X 
and p i — > Xp. If A is allowed to take on large values, there are only nine combinations of (9,p) which merit 
attention, namely all possible combinations of 9 = {0, 1, oo} and p = 0,p > 0,p < 0. These nine cases and 
their associated fixed points will also be the basis of analysis in subsequent sections of this paper. However, 
the cases corresponding to p < are physically less relevant, since eqn ( |4.1| ) then assigns the largest rates to 
the smallest cluster sizes. Figure p] shows schematically the effect of the contraction on the parameters 9, p. 



X + j) ~ X\oga + p \og(nX — X + x)dx 

Jo 



- 1 + fl-n) log 1 



(4.14) 
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VIII 
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Figure 1: The effect of the coarse-graining contraction or, equivalently, the renormalisation mapping, on the 
parameters 9,p in the Becker-Doring model with power-law rate coefficients. The nine fixed points of the 
map are denoted I-IX, and represent possible universality classes for the dynamics. 



In phase plane terminology, II has the form of an unstable node, while I, III, V, VIII are saddle points, 
since they act as attractors in one direction and repellors in the other (although they are at the limits of the 
allowable domain, so only have trajectories on one side of the critical point). Cases IV, VI, VII, IX are stable 
nodes since they act as attractors in both the 9 and the p directions. 
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5 Large-time asymptotics of the nine universality classes 



The previous section identified nine special, or limiting, cases which merit more detailed analysis, since 
following a coarse-grained contraction with large mesh parameter, any model will be approximated by one of 
these cases. To simplify the analysis we shall always choose a scaling such that x\ = 1, with a and (3 equal 
to zero or unity. In this section we present large time and large aggregation-number asymptotics. These are 
developed by means of a continuum limit, which replaces the system of ordinary differential equations by a 
single partial differential equation, and draws on the results derived by King & Wattis We apply these 
results in particular to our nine limiting cases. In a further paper these results will be generalised to systems 



whose rates are perturbed away from the simple polynomial expressions considered here [34 



5.1 Case I: pure fragmentation at a constant rate 

This corresponds to constant fragmentation with rate f3 n = 1 and no aggregation. Hence we expect the 
system to approach the equilibrium configuration x n = for all n > 2. The governing equation is 

x n = x n+ i - x n , for n > 2. (5.1) 

If the system is started from the initial conditions x n (0) = for n > 2, then the system does not change. All 
other initial conditions will approach this state as t — > oo. 



5.2 Case II: aggregation and fragmentation at constant rates 

This case corresponds to constant rate coefficients but with both coagulation and fragmentation present, and 
exactly balancing each other, a n = 1 = j3 n . Thus we expect the system to approach the equilibrium solution 
x r (t) = x\. The time evolution follows the equation 

x n = x n _i - 2x n + x n+ i, for n > 2, (5.2) 

which has the form of a "discrete diffusion" equation. 

In the limit of large time and large n, the discrete system Q5.2|) can be approximated by the diffusion 
equation 

I = JR. < 5 - 3 > 
at on 2 

whose solution is x n = erfc(n/2\/t). This equilibrium is approached by a purely diffusive mechanism, there 
being no advective term in ( |5.3| ), in contrast to Cases I and III (below) where advection is present. 

5.3 Case III: pure aggregation at a constant rate 

This case corresponds to constant aggregation with rate coefficients a n = 1 for all n, and no fragmentation 
(P n = for all n). Thus we expect the steady-state solution x n = x\ to be approached, which has steady-state 
flux L = ax n x\ = 1. The governing equations 

x n = x n _i - x n , for n > 2, (5.4) 

are exactly solved by 

n-2 ,k 

x^l-e-'Y: Tr (5-5) 
k=o K - 
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However it is not straightforward to see which mechanisms are driving the growth process when the solution 
is written in this form. An asymptotic approximation gives a much more intuitive idea of the dynamical 
behaviour. To describe the large time and large aggregation number asymptotics of this case, we approximate 
Qby 

dx x d 2 x dx .p. 
dt 2 dn 2 dn ' 

from which we can deduce the existence of a front of matter which propagates to larger aggregation numbers 
with speed unity, the front itself spreading out over a range of aggregation numbers proportional to t 1 ^ 2 . This 
can be shown by transforming to the new independent variable z = n — t yielding || = |fpr, which has the 
similarity solution x = ^erfc(z/ ^/t). Thus in the limit t — > oo 




x n = ^erfc [ —/=■ ) , (5.7) 

which is an asymptotic approximation of 



5.4 Case IV: pure fragmentation with rates increasing with cluster size 

In this case the only process occurring is fragmentation, and the fragmentation rate depends on cluster size, 
with larger clusters losing monomers at faster rates. So, as in Case I the system will tend to x n = for all 
n > 2. The system is governed by 

x n = n pX x n+1 -{n- l) pX x n , n>2. (5.8) 

Since the aggregation rates are greater than in Case I, the equilibrium solution will be approached more 
rapidly than in Case I. Matter is advected to smaller aggregation sizes, and removed from the system in 
monomeric form. 



5.5 Case V: aggregation and fragmentation with rates increasing with cluster 
size 

In this case we have both coagulation and fragmentation and both rates increase algebraically with aggrega- 
tion number. The system is governed by the equations 

x n = {n- l) pA x n -i - (n - l) pX x n - n pX x n + n pX x n+l . (5.9) 

The equilibrium solution is x n = 1 for all n, and by analogy with Case II (which also has 9 = 1) we might 
expect the large-time asymptotics to approach this solution. However, if we seek solutions from the wider 
class of steady-state solutions which are simply characterised by constant flux L (not necessarily equal to 
zero), then we find a family of solutions 

n— 1 i 

Xn = l-Lj2j^- (5-10) 

k=l 

From this family, we seek the solution which has the fastest possibly decay as n — ► oo. This occurs for 
L = l/((p\) (where ((z) = Y^T=i k~ z ^ s ^ e Riemann zeta function). Thus for p < 1/A, we have that L = 
and the equilibrium solution determines the large-time asymptotic behaviour of the system; however, for 
p > 1/A there is a steady-state solution (with L > 0) which has a faster decay at large r; hence the evolution 
will approach this solution rather than the equilibrium solution. 
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We shall now discuss the manner in which the equilibrium solution x n = 1 is approached in the case 
p < 1/A. The continuum limit for n ^> 1 is 

*=„P»f|£ + (5.11) 

at \on z n on J 

which immediately yields a similarity solution with variable n = n/t l '^ 2 ~ pX \ Assuming x = f(rj) with 
/(0) = 1 implies 

= /-^ A e X p(-^V(2-pA) 2 )^ 

/o 00 w-f A exp(-M 2 -fV(2-pA) 2 )^" 1 ' J 

In this case, there is no travelling diffusive wavefront; rather the balance of coagulation and fragmentation 
means that the equilibrium solution is reached by a purely diffusive mechanism. However, the similarity 
variable n shows that the aggregation number at a given concentration increases at a rate n oc t l ^ 2 ~ pX \ which 
is faster than the usual spreading rate of n oc y/t, this is due to the size-dependent coefficients accelerating 
mass transfer in the large n region. The precise details of the special case p — 1/A are not covered by this 
general solution, but the special case p = (Case II) can be obtained from the general theory for p < 1/A. 
When p > 1/A the analysis can be modified for the approach to the steady-state solution 



1 00 1 

^ SS = — -Et^t- (5.13) 



n C(PA) tn kPX 

The manner in which such a solution is approached can be found by writing x n (t) = x s ^ s tp n (t). We are then 
interested in the manner in which ip n — ► 1 as t —>■ oo. This is governed by the equation 

= f (n-l) pX + ^n-l - - (n-l)^n + U X - A] (5.14) 

\ X n J \ X n / 

where L is the steady-state flux (L = l/((p\)). For large n and at large times, it is appropriate to take the 
continuum limit, where x s ^ s can be approximated by l/(p\ — l)n pX ~ l ((p\), yielding 

1 d± | / 2~pA \ d± (5 15) 

n pX dt dn 2 \ n ) dn 

This equation possesses a similarity solution of the form ^(r,t) = f(n) where n = r/t 1 ^ 2 ^^ and 

u pX ~ 2 exp{-u 2 - pX /{2-p\) 2 )du 

/W = -> • (5-16) 

/ u pX ~ 2 exp{-u 2 - pX /(2-pX) 2 ) du 
Jo 

However this solution is only defined for 1/A < p < 2/A and is not valid for the case p = 1/A. In the special 
case p = 2/A equation ( |5.15|) is solved by 



Cleric (*^). (5.17) 



Although different states are approached in the cases p < 1/A and p > 1/A, in both cases the asymptotics 
are governed by a similarity solution of the form r] = n/t 1 ^ 2 ^ pX \ 
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5.6 Case VI: pure aggregation with rates increasing with cluster size 

This case corresponds to pure aggregation, the rate of aggregation growing with cluster size. There is thus 
no equilibrium solution; instead the system approaches the steady-state solution x n = n~ pX according to 

x n = (n- l) pA a; n _i - n pX x n n > 2. (5.18) 

The substitution ip(n,t) = n pX x n enables progress to be made. From the initial conditions ip n — for n>2, 
the steady-state solution ip n = 1 will be reached as time progresses. For p > 1/A the system is ill-posed due 
to the instantaneous transport of material to large aggregation numbers. For p < 1/A we assume that the 
equilibrium state is reached by a front moving from n — 1 to large n as time increases. To find the form of 
this wave at large times, we take the large n and large t limits, replacing eqn (|5.18j ) by 



^t- n P^(l^±^^±) Gig} 
dt [*dn* On)' lb ' i9j 

We define the position of the front by n = s(t), and transform to a moving coordinate frame z = n — s(t). 
We now show that the two 'outer' solutions ip = 1 valid for z — > — oo (that is n/s(t) < 1) and ip — which 



holds in the limit z — > +oo (that is n/s(t) > 1) can be joined smoothly by a transition layer. Equation (|5.19|) 
implies that ip(z,t) is determined by 

S-^Hr (*£-£)■ 

The leading order terms are those involving giving s = s pX , which is solved by 

s(t) ~ [(l-pA)t] 1/(1 " pA) , as t -> oo. (5.21) 

This result can be verified by considering the density, £ = Z^^Lx nx n ~ ^ 1_pA ~ (1 — s 2 ~ pX )/(p\ — 2), 
which implies that the density asymptotically increases according to g ~ s 1_pA s. However, the formula for the 
density also satisfies g = 1 + Z)^! n pX x n = s. Thus this calculation also gives s = s pX , confirming equation 

The first correction terms in (|5.20|) are 



1 dip 1 d 2 ip p\z dip 

g pX ~dt ~ Yd£ ~ ~~dz~' ^' ^ 

Since s(t) is monotonically increasing, we simplify further calculations by using s as our temporal variable. 
The next stage of the calculation is straightforward if p < 1/2 A, since in this case z ~ y/s and (|5.22|) reduces 
to |f = igf. Thus we find 



x n = \n~ pX ^ C ( in ~ S) ^~^ ) ■ (5-23) 

For p > 1/2 A, advection of the front dominates the diffusion processes and so the front retains its memory 
of the initial conditions. An alternative way of visualising this process is that the front moves so fast that it 
reaches n = oo before it has had time to adjust its shape to that of the erfc function. Thus we cannot specify 



in detail the shape of the front. More details of the special case p = 1/2A can be found in [37 
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5.7 Case VII: pure fragmentation with rates decreasing with cluster size 

Here there is no aggregation, so fragmentation will cause all clusters to break up, and as in Cases I and IV, 
the system will approach the state where all matter is in monomeric form (x n = for n > 2). Since the 
fragmentation decays rapidly with size, large clusters will take a very long time to break up. The determining 
equations are 

x n = n pX x n+1 - (n - l) pX x n , n>2. (5.24) 
Since fragmentation is weaker than in Cases I and IV, Case VII takes longer to reach equilibrium. 

5.8 Case VIII: aggregation and fragmentation with rates decreasing with clus- 
ter size 

This corresponds to coagulation and fragmentation being equally balanced, both decaying as the aggregation 
number increases. As in Cases II and V (with p < 1/A), the system approaches the equilibrium solution 
x n = 1 in the large-time limit. The similarity solution valid for p < 1/A derived in Case V is valid for this 
case too. The balance of coagulation and fragmentation forces means that the equilibrium solution is reached 
by a predominantly diffusive process. 



5.9 Case IX: pure aggregation with rates decreasing with cluster size 

The final case has no fragmentation, only an aggregation term which decays with increasing particle size. 
The system thus has no equilibrium solution, and will evolve to the divergent steady-state solution x n = n~ pX 
according to 

Xn = {n - l) pA x n _i - n pX x n , n>2. (5.25) 

The substitution ip n = n pX x n yields n~ pX ip n = ip n -i — VVi> an d progress can be made by forming a continuum 
approximation 

n -p\^P_ = 1 ,g 

dt 2 dn 2 dn 

As in previous pure aggregation cases, the large-time asymptotic approach to the steady-state is governed by 
a diffusive wave moving in the n-domain from n = 1 to large n; also as t increases the wavefront widens. The 
speed of propagation can be found by considering the density injected into the system, or by substituting 
n = s(t) + z into equation ( |5.26| ) and expanding in z. The leading order terms are 

1 d ^-^ + ^(s(tr px m-i- p —)- (5.27) 



s{t) pX dt 2 dz 2 dz 

At leading-order the dominant terms are those involving giving s = s pX since we assume z/s 1. The 
position of the front is given by 

s(t) ~ [(1 - p\)t} 1/{1 - pX) , as t — oo, (5.28) 

and we now use s as a new timescale to solve the next order terms in ([5727]), which yield ip = |erfc(z-\/((l — 
2p\)/(2s))). Thus 

x n (t) ~ ^- pA erfc \ [ — ~ 1 > as t ^ oo with n - s(t) = O (Vt) . (5.29) 
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From equation (|5.28|) it is seen that the effect of A is to reduce the exponent with which the position of 



the wavefront scales with time. Equation ( |5.29 ) shows that A increases the decay exponent of the steady- 
state solution from p to pX, and alters the scaling of the width of the wave. In the fully detailed system 
of equations the width asymptotes to \j2s/ (1 — 2p), which, under coarse-graining, one might expect to be 



mapped to y2s/(l — 2p)/A; however, solving the coarse-grained equations yields y2s/(l — 2pX) instead. 

6 Discussion 

We have described in detail the nine archetypal classes of behaviour into which the asymptotic dynamics of 
the Becker-Doring equations with power-law coefficients falls. These nine classes ( "Cases I-IX" of Section |5|) 
arise as a result of the application of a systematic renormalisation procedure to the fine-grained Becker-Doring 
equations and capture qualitatively different physical properties which are shared by all the models within 
the class, regardless of their microscopic differences. This amounts to a massive simplification of the original 
problem. In the present paper, we have concentrated on the detailed large-time asymptotic analysis of these 
nine cases. 

One specific detailed feature is worthy of note here. Our analysis of Cases IV, V and VI shows that 
coarse-graining retains the correct leading-order structure of the problem, though critical exponents of p = 1 
where the behaviour changes from fragmentation to aggregation dominated are mapped to p — 1/A. Thus if 
p < 1, to retain the correct qualitative behaviour one is limited to taking A < 1/p: taking the limit A — > oo 
causes the correct form of the kinetics to be lost, although the correct solution is still be approached in Cases 
IV and VI. In Case V, the equilibrium solution x n = 1 is approached if A < 1/p, and a steady-state solution 
if p > 1. However, if A > 1/p > 1, the coarse-grained system approaches a steady-state solution whereas the 
microscopic model approaches the equilibrium solution. 

In a companion paper we analyse the effect that random perturbations to the rate coefficients have 
on the above analysis. There we consider the modifications such perturbations make to the steady-state 
solutions, the kinetics of approach to these states and their effects on coarse-graining process itself. That 
analysis is necessary to confirm the existence and number of distinct universality classes in these Becker- 
Doring systems; as we have previously announced, all nine classes of behaviour may indeed be described as 
such (modulo, in Cases VII- IX, a particular requirement on the behaviour of the noise) [|16|1 . 

The renormalisation programme described here is of direct relevance to many of our recent publications on 
the application of generalised Becker-Doring models to a wide range of problems of physicochemical interest. 
A further paper currently in preparation [|35|] explores several of these nucleation and growth problems in 
which the renormalisation scheme is exact. 
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